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Abstract 

The present document aims at developing the formalism needed in order to describe 
the two-dimensional motion of a vehicle with Mecanum wheels, including the effects 
of resistive friction. The description of recently-acquired experimental data, on the 
basis of this model, is found satisfactory. 
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1 Introduction 

Despite their use in several applications, few efforts have been taken toward the 
theoretical description of the motion of Mecanum-wheeled systems. One recent 
article, dealing with the modelling of such systems [1], has not been edited 
to the extent of eliminating errors and misprints in the formulae it contains. 
A follow-up paper by the same authors [2] does constitute an improvement, 
yet it is not detailed enough in relation to the results obtained in that study. 
At the end their paper [2], the authors repeat a statement which appeared in 
an earlier version of the work [3], namely that 'More tests will be run on the 
single wheel tester to eliminate some of the errors seen in the current results'. 
However, the reader of these papers cannot (at least, in a direct manner) 
obtain a clear idea of the magnitude and seriousness of the 'errors seen in the 
current results' in terms of the description of the experimental data with the 
proposed model. 



The present paper serves two purposes: a) to provide a firm basis for the 
description of the motion of a Mecanum-wheeled vehicle and b) to establish an 
easy-to-implement solution, which could be used in the study of such systems. 

The mathematical basis for the description of mobile robots was developed 
by Muir and Neuman in the mid 1980s [4]; in that work, the authors had 
modelled the motion for six prototype wheeled mobile robots in terms of their 
so-called wheel equations. Despite its formality, however, the Muir-Neuman 
paper does not provide treatment of the effects of the resistive friction. 

The outline of the present paper is as follows. Section 2 will provide the details 
of the formalism which has been developed in order to model the motion of a 
Mecanum-wheeled vehicle on a horizontal plane; the motion will be determined 
in the general case, i.e., including resistive friction. Section 3 will confront the 
model of Section 2 with recently-acquired experimental data; the first part 
of the section will provide a short description of the experiments performed, 
whereas the second part will give the results of the analysis of the experimental 
data. Finally, the conclusions of the present work will be given in Section 4. It 
should be borne in mind, right from the beginning, that, as this study relies 
exclusively on very limited and specific data, it should not be regarded as 
complete; the proposed model must be confronted with experimental data, 
pertaining to arbitrary, not only translational motion. 



2 Method 

2.1 General 

The Mecanum wheel is frequently called 'the Hon wheel' after its inventor, 
Bengt Hon, who came up with the original design in 1973 while being em- 
ployed at the Swedish company Mecanum AB. It is a conventional wheel with 
a number of rollers attached to its circumference. The axis of rotation of the 
rollers is inclined with respect to the rotational 'plane' of the wheel, in a 
plane parallel to the axis of rotation (of the wheel). Omni-directional motion 
of Mecanum-wheeled vehicles is achieved by appropriately controlling the an- 
gular velocities of each wheel, as well as the direction of rotation (we will 
assume herein that each wheel is powered by a dedicated motor); in the gen- 
eral case, the motion is a combination of rotation and translation. If all wheels 
are turned in the same direction, at equal angular velocity, forward/backward 
motion of the vehicle is achieved. By rotating the wheels on the same side 
against each other, a sideways motion of the vehicle is achieved. There are 
velocity/rotation combinations resulting only in rotation of the vehicle. The 
mobile robot URANUS [5] (Fig. 1), the first flexible mobile robot designed 
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and constructed in order to provide a general-purpose mobile base to sup- 
port research in robot navigation, used Mecanum wheels for omni-directional 
motion. An overview of Mecanum-wheeled systems, albeit avoiding the math- 
ematical intricacies, including a discussion on their application, advantages, 
and disadvantages, may be found in Ref. [6]. 

2.2 Modelling of the motion of a vehicle with four Mecanum wheels: friction- 
less case 

Generally speaking, it does not sound logical to talk about motion of vehicles 
in the case of vanishing friction; it is through friction that motion is enabled. 
What is actually implied when using the term 'frictionless' in the description 
of a motion is that the tractive forces are large enough to enable motion, 
whereas the resistive ones are conveniently small so that they may safely be 
ignored. In this section, the formalism needed in order to describe such an 
ideal situation will be developed. 

Figure 2 provides a top view of a (rectangular) vehicle featuring four Mecanum 
wheels, along with its attached coordinate system (x,y), the origin of which is 
assumed to be the geometrical centre of the rectangle; the wheels are identified 
by the numbers 1 . . . 4, starting from the right-bottom corner (i.e., from the 
right-rear wheel of the vehicle) and proceeding in the counter-clockwise direc- 
tion. The angular velocities W1...4 are defined positive for translational motion 
in the forward direction (increasing y). 

The driving (motor) force (thrust) Fj acting on wheel % of the vehicle (chosen 
to be wheel 2), along with its decomposition into one force (F itP ) parallel to 
the rotational axis of the roller (which is in contact with the ground at that 
moment) and one in the transverse direction (F ijt ), are shown in Fig. 3. The 
angle between the transverse direction and the rotational plane of the wheel 
is denoted as a G [0,7r). (The quantity sin a is also known as the 'efficiency 
of the wheel'.) Since the rollers rotate freely around their axle, there is no 

— * 

traction along the transverse direction; therefore, the force F i)t can safely be 
ignored when studying the motion of the vehicle. The relation between F ijP 
and Fi (indicating the corresponding moduli of the two vectors) reads as: 
F itP = Fi sin a. 

Finally, the only relevant force, F i;P , may be decomposed into forces along the 
axes of the attached coordinate system (see Fig. 4). The geometry dictates 
that F i>x = F itP cos a = Fi sin a cos a and F ijV = F itP sin a = Fi sin 2 a. 

To obtain the expression for the total force acting on the vehicle in the x 
and y directions, one must consider the different orientations of the axles of 
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the rollers on each wheel; in the configuration which is 'casually' used, these 
orientations are identical for diagonal wheels and opposite (mirror images) 
for the two front/back wheels. From now on, the sign of Fj will be omitted 
(i.e., the tractive force will always be assumed positive) and the direction of 
applied force will instead be taken into account in the rotational direction of 
each wheel Q Assuming that the quantity a represents the orientation of the 
roller axes of wheels 2 and 4, one may easily deduce that the force acting on 
the vehicle in the x direction is: 

4 4 

F x = ^l F i,x = sin a cosa^(-l)*s^n(a; i )Fj ; (1) 

i=l i=l 

the force in the y direction reads as: 

4 4 

F y = F i,y = sin 2 a ^ sgn(ui)Fi . (2) 

i=l i=l 

It must be mentioned that Eqs. (9, 10) of Ref. [1] contain one additional 
quantity (Ki), introduced as 'the wheel constant, dependent on the number 
of rollers per wheel and [on] how tight the rollers are [fixed] on the wheel's 
hub'. Given that the goal herein is to develop the formalism for similar wheels 
(material and number of rollers, as well as inclination a), the K^s can be 
assumed constant and, as such, be absorbed in a redefinition of the tractive 
forces Fi. Finally, the total force acting on the vehicle is given by 

F = F X U X + FyUy , (3) 

where u x and u y denote the unit vectors along the x and y directions, respec- 
tively. In the special case that u\ = U2 = 0J3 = U4, the lateral forces cancel 
out (i.e., the force in the x direction vanishes, see Fig. 5). 

The motion of a rigid body in two dimensions is a superposition of a) a trans- 
lational motion of its centre of mass (CM) and b) a rotation around the axis 
passing through the body's CM, normal to the plane on which the motion is 
confined. 

To determine the translational motion of the vehicle in the laboratory, one 
needs to obtain the components of the force in that reference frame from 
those of the CM coordinate system. If the latter is rotated by an angle <fi with 
respect to the former (a parallel translation leaves the orientation angle of a 



1 The projection of the angular velocity along the x direction (simply denoted 
as uji) is either positive or negative. 
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vector invariant), the transformation reads as 

cos 4> — sin <j> . . 

. . ( 4 ) 

sin cos 

The total force Fl in the laboratory coordinate system is given by 

Fl = F x u x + F Y u Y , (5) 

where ux and uy denote the unit vectors along the X and Y directions. 

The translational motion of the vehicle in the laboratory coordinate system is 
obtained on the basis of Fl of Eq. (5) and of the mass m (also containing the 
equivalent mass of loads) of the vehicle, whereas the rotation will involve the 

— * 

torque T acting on the vehicle. 



f = X>1 x F t = X)(r1 -r^) x F t = £>1 x F t x ]T = 

i=l i=l i=l i=l 

4 

^rlxF.-r^xF, (6) 



i=l 



where r C u identifies the vehicle's CM in the attached coordinate system (x,y), 

— * 

fl the position vector of wheel % in (x,y), and r[ the position vector of wheel 

— * 

i in the coordinate system (x',y') (see Fig. 6). The torque T will generate a 
rotation around the z'-axis of the coordinate system, passing through the CM 
of the vehicle; the rotation obeys the equation 

f=m x , , (7) 

where / is the moment of inertia of the vehicle (around z'), <fi represents 
the angular acceleration, and u z i is the unit vector along the z' direction 
(u z > = u x i x u y /). 



2. 3 Determination of the position of the vehicle 



The problem of the determination of the position of the vehicle as a function 

of time t will be split up into two parts, a) Evaluation of the position of 

the vehicle's CM as a function of time (i.e., determination of the function 
> > 

RcM(t), where Rcm denotes the position vector of the CM in the laboratory 
coordinate system), b) Evaluation of the orientation angle of the vehicle (i.e., 
of the amount of the rotation around the z'-axis) as a function of time (i.e., 
determination of the function <f>(t)). A schematic form of the relation between 
these two coordinate systems is shown in Fig. 7. 
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Let us assume that the position of the CM of the vehicle at t — t\ is known 
(i.e., that Rcuiti) is given), as is the velocity Vc-m(^i)- An estimation of 
RcM(ti + At) and VcM(ti + At) will first be obtained. The time interval At 
is chosen small enough to ensure that the acceleration between t\ and t\ + At 
be safely considered constant. If A C m = F L /m, with F L taken from Eq. (5), 

one obtains 

Vc2(t x + At) = Vcjif (ti) + A CAf At (8) 

and 

i^J(t! + At) = R^rih) + Vc2(ti)At + -A^(At) 2 . (9) 

The corresponding expressions for the angular velocity and angle may easily 
be obtained by using the standard substitutions: F — » T, v — > Co, r — > <fi, and 
m — > I. Therefore, 

f 

w(ti + At) = w(ti) + yAt (10) 

and 

— * 

${h + At) = + ^(tx)At + (At) 2 . (11) 

2.4 Modelling of the motion of a vehicle with four Mecanum wheels: intro- 
duction of resistive friction 

Resistive friction is generally categorised as static and dynamic. 

• Static resistive friction results in the force between two objects in contact 
when they are not in motion relative to each other. The static friction 
prevents the start-up of motion, as long as the tractive force F m remains 
smaller than a fraction of the (modulus of the) normal force N exerted 
between the two objects, i.e., as long as F m < F a = fi a N; in this relation, 
ji a is identified as the coefficient of static friction. 

• Dynamic resistive friction occurs when two objects in contact are moving 
relative to (i.e., rubbing against) each other; the coefficient of the dynamic 
friction will be denoted as ji K . 

The coefficient /i K is usually considered to be smaller than \i a (for the same 
materials), though this condition does not qualify as a necessary theoretical 
constraint; throughout the present work, it will be assumed that /i a = /i K . 

One must distinguish between two types of resistive friction, namely sliding 
and rolling; the former pertains to translational motion, the latter to rota- 
tional. In both cases, no distinction will be made between static and dynamic 
effects. The coefficient of sliding friction will be denoted as fj, s , the one of 
rolling friction as /i r . 
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In addition to the aforementioned types of resistive forces, the viscous resistive 
friction (introduced in analogy to the force caused by the viscosity of lubri- 
cants) will be taken into account; it depends linearly on the relative velocity 
v between the two objects in contact: F v = —fi v Nv, where \i v is the viscous- 
friction coefficient. No distinction will be made between sliding and rolling 
viscous resistive friction. As a matter of fact, the introduction of the viscous 
resistive friction in the present work is mathematically equivalent to assuming 
that the sliding- and rolling-friction coefficients are not constants, but depend 
(linearly, and in the same manner) on the velocity v. 

Let us denote the velocity of wheel i in the laboratory coordinate system at a 
time instant t as Vi. Evidently, 



Vi =V CM + R(0) 



u x r' : 



Vcm - <j> (y'i cos + x[ sin 4>)ux + (y'i sin <j> — x[ cos 0)«y , (12) 

where R(0) denotes the 2x2 rotation matrix appearing on the right-hand 
side of Eq. (4). 

— * — * 

Each vector Vi may be decomposed as: Vi = V^xUx + ViyUy. To determine 
the direction of the resistive force on wheel i, one must take three factors into 
account: a) the orientation of the vehicle in the laboratory coordinate system 
(angle 0), b) the orientation of the axle of the rollers with respect to the 
rotational plane of the wheel, and c) the direction of rotation imposed on the 
wheel. For wheels 1 and 3, the orientation angle (of the axes of the rollers) in 
the laboratory coordinate system (i.e., the angle with the X-axis) is equal to 
7r — a + 0; for wheels 2 and 4, it is equal to a + 0. Therefore, the directional 
vector for wheels 2 and 4 is Ui tP = cos(o; + <p)ux + sin(o; + 0)«y; for wheels 1 
and 3, u^ p = cos(7r — a + <p)ux + sm(ii — a + 0){ty. The projection of the vector 

— * 

Vi along the appropriate direction will determine the sign of the components 
of the resistive force (we will deal with its modulus | fa \ shortly), namely of 

— * 

fi t x and fi t y- For wheels 2 and 4, the projection of the vector Vi along the Ui tP 
direction is equal to V^ p = V^x cos(a + 0) + V^y sin(o; + 0); for wheels 1 and 
3, Vi tP = Vi t x cos(7r — a + 0) + Vi t y sin(n — a + 0). Therefore, if a G (0, |), 
the two components of the resistive force have the same sign (i.e., either both 
positive or both negative) for wheels 2 and 4; for wheels 1 and 3, they are of 
opposite sign. 

For the evaluation of the modulus of the resistive force, one must first deter- 
mine the corresponding normal force iVj acting on wheel %. For planar motion 
of the vehicle, 

Ni = N L *~ \< \ L v-\ti\ ? (13) 

L X Ly 

where L x and L y are the distances between the wheels in the x and y directions 
('track' and 'wheelbase' of the vehicle, respectively); the coordinates of wheel 
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i (x'i and y[) correspond to the CM coordinate system (see Fig. 6). If rem = 0, 
all forces iVj have the same magnitude (i.e., one quarter of the total load N). To 
obtain /j from iVj, one must multiply by the appropriate coefficient of friction, 
i.e., either /x s (if the wheel % is not moving in the direction of the traction) or 
\i r (if it is); if the wheel is in motion, the viscous- friction force must also be 
included. Finally, the forces F x and F y of Eqs. (1, 2) will be redefined after 
taking into account the contributions of the resistive forces. 

At this point, one remark is due. Rewriting Eq. (1) and absorbing each 
sgn(ui) factor in the corresponding Fi, one obtains: F x = (— F\ + F 2 — F 3 + 
F 4 ) sin a cos a; Equation (2) leads to: F y = (F x + F 2 + F 3 + F 4 ) sin 2 a. Assum- 
ing no offset of the CM of the vehicle with respect to the attached coordinate 
system (x,y), there is no rotation of the vehicle if the resulting net torque 
vanishes; it may easily be shown that the condition for no rotation is: 

sin a(F 1 + F 2 — F 3 — F 4 )(L X sin a — L y cos a) = . (14) 

This equation admits three solutions, one of which (namely, sin a = 0) is trivial 
and uninteresting from the point of view of a practical application. However, 
the two remaining conditions are worth noting. Evidently, in case of right/left 
balance of the forces (i.e., when Fi + F 2 = F 3 + F 4 ), the net torque (in the 
standard configuration of the Mecanum wheels) vanishes. The last condition 
is the most interesting one: if the track and the wheelbase are chosen in such 
a way as to satisfy the condition tana = L y /L x , the resulting net torque 
always vanishes irrespective of the magnitudes of the applied forces and of the 
direction of the rotation of the Mecanum wheelsl (Given that both L x and 
L y are positive, the last equation becomes relevant only when a G (0, 7r/2).) 
This result must be borne in mind in case that the rotational effects must be 
suppressed. Finally, one may show that, if the condition Fi + F 2 = F 3 + F 4 
is fulfilled, the vehicle will move along a straight path, the direction of which 
will involve the angle tan" 1 ( ^^-f^ ~ ) + 0(0) with respect to the X axis of 
the laboratory coordinate system. 

2. 5 Derivation of the solution 

Let us assume that, at time ti, the known quantities are: 

• the position vector of the CM of the vehicle Rcm, 

• the velocity of the CM of the vehicle Vcm, 

• the orientation angle 0, 

• the angular velocity u, and 

• the tractive forces F { for each wheel. 

The output at time t\ + At comprises: 
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• the position vector of the CM of the vehicle, 

• the velocity of the CM of the vehicle, 

• the orientation angle, and 

• the angular velocity. 

Compared to systems using standard wheels, it is more complex to include 
the resistive effects in the general motion of a Mecanum-wheeled vehicle. In 
the former case, the wheels move in a concerted manner, thus creating either 
forward or backward overall motion; a vehicle using Mecanum wheels is bound 
to behave differently. Not only can the wheels be set to rotate at different an- 
gular velocities, but they might also be set to rotate against each other. As a 
result, the general state of motion of each wheel is best described as 'rolling 
with sliding'; this implies that, in the general case, the sliding- friction coeffi- 
cient appears to be the appropriate one to use in the expressions modelling the 
motion. On the other hand, if the motion is concerted (i.e., when the wheels 
are set to move in the same direction, at equal velocities), the coefficient of 
rolling resistive friction should better be used. To enable the inclusion in the 
description of the motion of the effects pertaining to friction and simultane- 
ously retain simplicity, the following scheme is proposed. 

(1) Loop over the four wheels: 

• Estimate Vi via Eq. (12) and its projection along the u^ p direction, 
i.e., the V^ p component entering the evaluation of the resistive friction 
exerted on the wheel by the floor. 

• Estimate the component of the thrust along the Ui tP direction (i.e., F itP ). 

• If | Vi iP | exceeds a tolerance limit (denoted in the flowchart below by the 
quantity e), then the resistive force is given by the sum of the rolling- 
and viscous-friction contributions (fi iP = —sgn(Vi jP )ii r NiU itP — ii v NiVi tP ) 
if Vi tP and F itP are of the same sign, whereas by the sum of the sliding- 

— * — * 

and viscous-friction contributions {f^ p = —sgn(Vi :P )[i s NiUi )P — [i v NiVi tP ) 
if Vi tP and Fj iP are of opposite sign; if, on the other hand, | Vi jP \< e 
but | F i>p |> fi r Ni, then the friction is overcome by the thrust and the 
condition for the motion of the wheel along the Ui jP direction is fulfilled 
(the motion will be allowed within the current period of time At). Of 
course, if the thrust has not overcome the friction, the wheel will not 
move within the time interval At. 

• Redefine F itP by taking the resistive forces f^ p into account. The follow- 
ing flowchart describes this simple procedure, applied to each wheel % 
of the vehicle. 

— * 

(2) Evaluate the components F x and F y from the forces F i>p (by consider- 
ing the projections of the unit vectors Ui iP onto the x and y directions); 
evaluate Fx and Fy from the matrix equation (4). 

(3) Use expressions (8-11) to obtain the output values of the relevant physical 
quantities at time t\ + At. 
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TP. TP. — 

± i,p ~ ± i,p 

sgn(F itP )n r NiU i)P 



yes 




yes 



W. ^ TP. - 

- 1 i,p ^ - 1 i,p 

sgn{Vi,p)yL r NiU ijP 



It is convenient to place the CM of the vehicle at t — at the origin of the 

> — % 

laboratory coordinate system, at rest; therefore, Rcm(0) = and Vca/(0) = 0. 

Additionally, 0(0) = and cJ(0) = 0. One may then start to apply the tractive 

forces to the wheels, perhaps using a linear model for t < T, and the 

nominal values Fi for t > T; a reasonable value for T may be a few tenths of 

Is. Concerning the time increment At, it must be sufficiently small to ensure 

the accuracy of the output. Assuming that the execution time is not an issue, 

one might put forth a scheme of subdividing the original time increment At 

into segments 5t, obtain the relevant physical quantities (at ti + At) for several 

St values, fit the values with a (simple) polynomial, and (finally) obtain the 

results at t\ + At using the extrapolated values to St = 0. 
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2. 6 Simple examples of motion 



The implementation of the formalism, developed in Sections 2.2-2.5, led to the 
results of Table 1 for a few simple motions; a similar table may be found in 
Ref. [1] 

A change of orientation of the vehicle y axis (with respect to the laboratory 
coordinate system) may be caused not only via the application of different 
thrusts (or rotations) on the Mecanum wheels, but also by shifting the CM of 
the vehicle (away from its geometrical centre). This is evident after inspecting 
Eq. (6); if tcm ^ 0, a (non-zero) torque is generated as soon as the vectors 
t C m and F are not aligned (or anti-aligned). The results of Table 1 have been 
obtained with r C M = 0. 

One case of motion of a 16.6kg vehicle is displayed in Figs. 8-10, over a period 
of 4s, in 100ms steps. The parameters generating this motion are (see next 
section): v\ = 10, v% = 12, u 3 = 14, z/ 4 = 16 revolutions per minute; sgn{u\) = 
sgn(u2) = sgn(cjs) = sgn{u)i) = +1 ('++++' in the notation of Table 1); 
N = 250kp; L x = 40, L y = 58cm; vqm = (5u x + 5u y )mm; a = 130°; fi s = 0.51, 
fi r = 0.084; fi v = 0.15s/m. The starting position in the laboratory coordinate 
system (X, Y, (p) is (20cm, 25cm, 10°). The thrusts are applied to the wheels 
linearly over T = 300ms; the nominal thrust values are assumed to have been 
reached at t — T (see last paragraph of the previous section). 

In the present work, the assumption is that the load is applied at one spe- 
cific point of the vehicle surface. As a result, the solution (extracted for the 
rotation of the vehicle) is, generally speaking, not very realistic; in reality, the 
application of the load involves the contact of two surfaces. Consequently, the 
amount the vehicle rotates in the ideal case is expected to exceed the actual 
amounts obtained in realistic cases. This is due to two reasons. Firstly, if the 
rigid object (which is hard to move in the attached coordinate system (x,y)) is 
placed onto the Mecanum-wheeled vehicle, part of the driving torque will be 
'wasted' in balancing the frictional torque, in analogy to the trade-off between 
tractive and resistive forces in translational motion. Secondly, the value of the 
moment of inertia I in Eqs. (7, 10, 11) will have to include an unknown (pos- 
itive) contribution pertaining to the geometrical characteristics of the object 
comprising the load. One way to solve part of the former problem is to allow 
the rotation of the vehicle only if the driving torque exceeds a specific value, 
and subtract a (perhaps, different) contribution from the driving torque after 
the rotation of the vehicle has been initiated (i.e., distinguish between static 
and dynamic frictional torques); without specific information on the charac- 



2 One must bear in mind that Ref. [1] differs from this work in some definitions, 
e.g., in the sign convention of the Wj's. 
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teristics of the object responsible for the load, it is not possible to obtain a 
solution to the latter problem. 



3 Validation of the model 

3. 1 Experiments 

A number of tests have been carried out in order to validate the model of 
Section 2 and to examine whether the values of its parameters could be de- 
termined from experimental data; the last of these tests (performed on the 
concrete floor of the IMS laboratory on October 3, 2012) was documented. 
The results of the previous tests had already shown that: 

• the vehicle was controllable, in fact, equally- well controllable over straight 
paths set at angles between and 90° with respect to the X axis of the 
laboratory coordinate system (initial condition: 0(0) = 0°), and 

• the angular uncertainty in its motion did not exceed 1° (root-mean-square 
value of the difference between expected and actual directions). 

For the documented data taking, a number of angular-velocity combinations 
(as well as rotational directions), resulting in a motion of the 16.6kg prototype 
within the first quarter-circle (i.e., X(t) > and Y(t) > 0), were applied to 
the wheels of the vehicle; the coordinates of the end points (of the trajectories) 
were measured relative to the walls of the test room, using a laser distance me- 
ter (Leika DistoTM A2 Laser). The motion was interrupted when the distance 
between the starting and final positions of the right-rear corner of the vehicle 
was about 2m (due to imperfection of the floor, the motion was interrupted 
earlier, at distances between 1.6 and 1.7m, below 20°). The details on these 
data may be found in Table 2. (The inclination angle a in the prototype is 
135°.) 

The friction coefficients were also measured in the two following cases: a) the 
brakes were applied unto the wheels and b) the wheels were allowed to rotate 
freely. The former measurement yields the coefficient of the sliding friction, 
whereas the latter one fixes the coefficient of the rolling friction. The results 
were: fi, = 0.51 ± 0.03 and (i r = 0.084 ± 0.006. 

The experimental data (a total of M = 32 measurements of the coordinates of 
the end points of straight paths, as well as an equal amount of data pertaining 
to the time needed in order to accomplish each motion) were submitted for 
optimisation. Two additional parameters were introduced (denoted as A and 
B), to associate the angular velocity Ui of each wheel with the tractive force F { 
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acting on that wheel P~[ a simple linear model was put forth: Fi = A Ru>i + B , 
where R stands for the radius of the wheels used in the prototype (R = 5cm). 



For a given set of parameter values, the motion of the vehicle was deter- 
mined until the measured end point was 'boxed'. An estimation of the short- 
est distance d{ of the resulting trajectory to the measured data point was 
next obtained. The uncertainty in this determination Sdi, to be assigned to 
each experimental value, was assumed proportional to the distance between 
the starting and final positions (of the reference point of the vehicle); the 
assumed angular uncertainty of 1° is translated into a spatial uncertainty of 
about 3.5cm over the distance of 2m. As the distances between the starting 
and final positions were kept close to 2m, the assigned statistical uncertainties 
of the input data (above 20°) were almost equal for the entire set of measure- 
ments. The optimisation was achieved via the minimisation of the x 2 function, 
defined by the formula: 



times corresponding to the particular motion; the uncertainties 8ti were 'glob- 
ally' set to 2s (which may have been an overestimation of the true uncer- 
tainties). The free parameters were varied until the description of the input 
experimental data was optimised (i.e., until the \ 2 minimum was obtained). 
Given that two physical quantities (i.e., the distance di and the time t1 xp ) per 
case were optimised, the total number of input data points is 2M = 64. 

For the purpose of the optimisation, the standard MINUIT package [7] of the 
CERN library was used, namely the C++ version of the library [8]. Extensive 
information on this software package may be obtained online [9]. The latest 
release of the code is available from Ref. [10]. Version 5.28.00 of Minuit2 was 
incorporated in the analysis framework. Each optimisation was achieved on 
the basis of the (robust) SIMPLEX-MINIMIZE-MIGRAD-MINOS chain. 

The first fits to the experimental data, treating all five model parameters 
(i.e., the three coefficients of friction, A, and B) as free, revealed the presence 
of strong correlations. To mitigate these correlations, a decision had to be 
made as to which parameters could be fixed from external sources. Given that 
the coefficients of the sliding and rolling friction had already been obtained 
from direct measurements, it was decided to use these two results and vary 

3 In order to determine the motion of the vehicle, it is necessary to associate the 
angular velocity uii of each wheel with the tractive force Fi acting on that wheel. In 
the future, this relation must be established after analysing dedicated experimental 
data, corresponding to the motors, the materials involved, and the geometrical and 
physical characteristics of the specific wheels used in the vehicle. 




(15) 



where tf p and tf 1 denote the measured and fitted (estimated with the model) 
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only fi v , A, and B in the fits. Due to the largeness of the correlations, a 
cautious stepwise approach during the next steps of the optimisation phase 
was followed. Apart from occasional fixing and releasing some of the model 
parameters (thus enabling the examination of the sensitivity of the results 
to the variation of the model parameters, as well as the assessment of the 
performance of the optimisation algorithms in the problem), the optimisation 
scheme was also run including and excluding the measurements relating to the 
duration of the motion; the results of this test demonstrated that the problems 
encountered in the data analysis did not originate from the inclusion of the 
time measurements in the fit. However, despite all these efforts, no satisfactory 
results were obtained; a number of fits failed to yield meaningful output (e.g., 
the Hessian matrix could not be evaluated) as the optimisation algorithms 
were trapped in spurious y 2 minima. (A dependence of the results of the fit 
on the initial guess of the model-parameter values was also observed.) On 
almost all occasions, the free (unconstrained) fits drifted as a result of the 
existing correlations. In order to cope with this behaviour, it was subsequently 
decided to use available information on some technical details of the hardware 
being used in the tests, and impose appropriate constraints on the model 
parameters (e.g., pertaining to the maximal power and torque of the motors 
used in the prototype), on top of the obvious physically-imposed constraints 
(e.g., positivity of the friction coefficients, etc.). Despite these actions, the 
tendency of the fits to drift persisted. 

It currently seems that the only way to obtain a stable solution was to perform 
the fits at fixed \i v values, taken from a representative interval. The question 
of what 'representative' means is better answered after considering the viscous 
friction as the velocity-dependent component of the (varying with the velocity) 
sliding/rolling friction force [^j From this viewpoint, one could set limits for 
the allowed variation of the sliding-friction coefficient between two reference 
velocities, say, between and lm/s (arbitrary choice). It should suffice to 
assume a variation of 100% (of the sliding-friction coefficient, assumed to be 
velocity-dependent) between these two reference velocities. In that case, the 
100% variation of the coefficient is equivalent to varying the parameter /i„ 
from to the value of fi s (in s/m units). Since each optimisation (at a fixed 
fi v value) involves two free parameters, the number of degrees of freedom in 
the fits was NDF = 64 - 2 = 62. 



4 In other words, whichever the dependence of the sliding/rolling friction on the 
relative velocity of the two surfaces in contact, it will be modelled via a linear func- 
tion; given the smallness of the typical velocities in the problem, this approximation 
is expected to be good. 
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3.2 Results 



The results of the optimisation for fixed values of the parameter /i v from (0.05 
to 0.50s/m, with a step of 0.05s/m) are shown in Table 3. Inspection of the 
table (and after trivial numerical operations on the data it contains) provides 
an answer to the question of the failure of the optimisation scheme in yield- 
ing meaningful results. The y 2 function appears to be monotonic with \i v (at 
least for the fi v interval of interest, as chosen herein). This behaviour of the 
minimisation function prevents the optimisation algorithms from terminating 
successfully. The fit drifts in search of the x 2 minimum, which (judging from 
the corresponding plot) is nowhere in the vicinity of the interval of represen- 
tative \i v values. In any case, Table 3 will surely be useful when simulating 
the motion of the vehicle in the general case; the user may choose any of the 
shown (/i v ,A,B) combinations, exempting perhaps only the first row (i.e., the 
record corresponding to /i v = 0.05s/m). 

In view of these results, it makes sense to restrict the analysis to a reasonable 
interval of \i v values, e.g., 0.10 < /i v < 0.20s/m. In that case, the average of 
the three minimal x 2 values is about 105.9 for 62 degrees of freedom in the fit, 
i.e., about 1.7 units (on average) per degree of freedom; this reduced-^ 2 value 
is reasonable. A more careful analysis of the results of the fits indicated that 
the motion was more controllable for small (say, below about 35°) and large 
(above about 75°) angles; the data in these regions were better reproduced 
(than those outside) with the model of Section 2. The worst cases in the 
reproduction of the coordinates of the end points occurred around 60°, where 
the reduced-^ 2 values were as large as 6.7, showing a deterioration from large 
to small angular velocity. This might be due to the influence of effects which 
have not been taken into account in the model of Section 2 (e.g., see Ref. [2]). 

The reproduction of the time duration of the motion (see Fig. 11) was found 
superior to that of the coordinates of the end points of the trajectories. The x 2 
value, corresponding to the reproduction of the 'temporal' data, was equal to 
only 28.0 units for 32 measurements. The linear-correlation coefficient between 
measured and fitted values exceeded 0.9975. 



3.3 Further comments 

It is true that another approach could have been followed in the optimisation 
phase. For instance, the trajectory could be determined for a time equal to 
the one measured in the experiment; the coordinates of the end points of 
the resulting trajectory could then be compared directly to those obtained 
experimentally. In this scheme, only the first term of the minimisation function 
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of Eq. (15) is relevant. One could also introduce the time uncertainty Sti, by 
appropriately translating it into a spatial uncertainty and combining it (e.g., 
quadratically) with the (assigned) spatial uncertainties 5di of the data points. 
In principle, this approach and the one followed herein should lead to similar 
results. 

The measurements obtained during the experiment rely on the determination 
of the distances of the reference point of the prototype to the walls of the test 
room. Concerning these measurements, the assumption is that the planes of 
these walls are perpendicular to each other (Cartesian laboratory coordinate 
system). Unfortunately, a deviation of about 1° was observed after simply 
placing a right-angle ruler at the corner of the two walls. Therefore, it makes 
sense to take the numbers of Table 2 with a grain of salt. Equivalently, one 
could accept these values (and the bias they might introduce into the analysis), 
but augment the input uncertainties of the data points; the latter strategy 
was adopted herein. (The use of a laser scanner, as in Ref. [2], would be a 
considerable improvement at this point.) 

The model of Section 2 was put forth in order to understand the motion 
of one vehicle using Mecanum wheels, so that a more complex system (e.g., a 
swarm of such vehicles) be efficiently studied. In this context, it is not relevant 
whether all physical effects have analytically been included in the proposed 
model, as long as the influence of any missing pieces is captured by the model 
parameters, which (in this case) lose their clear physical interpretation and 
become effective. Finally, by no means should the tests performed herein be 
considered complete; to thoroughly test the model of Section 2, one would 
need to develop a test platform similar to the one shown in Fig. 7 of Ref. [2] . 



4 Conclusions 

In the present paper, given are the details of the modelling of the omni- 
directional, two-dimensional motion of a vehicle using Mecanum wheels. The 
model is based on simple Mechanics and includes the effects of rolling, slid- 
ing, and viscous friction. The motion is obtained in the general case, i.e., for 
arbitrary thrusts applied to the wheels of the vehicle. 

The three friction coefficients comprise the parameters of the model. Two ad- 
ditional parameters were introduced at the present stage, in order to account 
for the dependence of the tractive force (applied to each wheel) on the cor- 
responding angular velocity; in the future, these parameters could either be 
fixed from dedicated experimental data or be directly derived from the de- 
tailed technical specifications of the motors, as well as the geometrical and 
physical properties of the wheels. 
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The model was confronted with recently-acquired experimental data, compris- 
ing the coordinates of the end positions of one reference point of the vehicle 
for straight paths in the first quarter-circle, as well as the times needed in or- 
der to achieve each motion. Given the strong correlations, which were present 
among the model parameters in the optimisation of the description of the 
experimental data, two of the friction coefficients (i.e., of rolling and sliding 
friction) were fixed from direct measurements. This action, however, did not 
suppress the correlations. To obtain meaningful results from the experimen- 
tal data, the fits had to be performed at fixed values of the viscous-friction 
coefficient, taken from an interval which appears to be reasonably broad. 

The present work may serve as the basis for the simulation of the motion of one 
vehicle (for arbitrary thrusts, inclination angle of the axles of the rollers, loads, 
and offsets of the centre of mass of the vehicle), as well as of the concerted 
motion of a swarm of such vehicles. Obviously, further work is needed both 
on the theoretical, as well as on the experimental aspects in this topic. 
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Table 1 



A few cases of simple motion of a Mecanum-wheeled vehicle; a similar table may 
be found in Rcf. [1]. CW denotes 'clockwise', CCW 'counter-clockwise'. The signs 
correspond to the quantities sgniuji); indicates no rotation. The cardinal directions 
correspond to a top view of the motion of the vehicle, i.e., E is the direction of 
positive X (at constant 1"). The same thrust value was applied to the wheels. The 
orientation of the rollers axles was set to 45°. 



Wheel 1 


Wheel 2 


Wheel 3 


Wheel 4 


Direction of motion 


+ 


+ 


+ 


+ 


N 





+ 


+ 





N 


+ 








+ 


N 


— 






— 




— 







S 
S 


— 








— 


s 


+ 





+ 





w 


+ 











w 








+ 





w 




+ 




+ 


E 




+ 








E 










+ 


E 


+ 





+ 





NW 












SE 





+ 





+ 


NE 












SW 


+ 


+ 






CW Rotation 






+ 


+ 


CCW Rotation 



19 



Table 2 



Details on the input data acquired on the concrete floor of the IMS laboratory on 
October 3, 2012. These values are analysed in Section 3. The quantity v denotes 
the revolution rate (number of revolutions of the relevant wheel per unit time); in 
all cases, 1/3 = v\ and 1/4 = v 2. 



tr p (s) 


Xr (m) 


Yr (m) 


v\ (min 1 ) 


v> 2 (min -1 ) 


26.6 


1.5770 


0.2615 


14.55 


-10.19 


35.2 


1.6180 


0.2675 


11.25 


-7.88 


49.7 


1.6190 


0.2755 


7.95 


-5.57 


85.0 


1.6320 


0.2665 


4.65 


-3.26 


31.7 


1.6130 


0.5855 


14.55 


-6.78 


41.1 


1.6440 


0.5805 


11.40 


-5.32 


56.8 


1.6500 


0.5755 


8.25 


-3.85 


91.1 


1.6380 


0.5785 


5.10 


-2.38 


36.8 


1.6270 


0.9525 


14.55 


-3.90 


47.0 


1.6500 


0.9715 


11.55 


-3.09 


62.0 


1.6440 


0.9585 


8.70 


-2.33 


92.6 


1.6550 


0.9575 


5.85 


-1.57 


50.1 


1.8820 


1.6515 


14.55 


-1.27 


55.0 


1.7440 


1.5405 


12.30 


-1.08 


69.9 


1.8430 


1.6165 


10.20 


-0.89 


86.5 


1.7970 


1.6025 


8.10 


-0.71 


46.5 


1.4120 


1.8425 


14.55 


1.27 


60.3 


1.4160 


1.8515 


11.25 


0.98 


83.5 


1.4400 


1.8915 


8.25 


0.72 


123.4 


1.5400 


2.0485 


6.00 


0.53 


44.9 


1.0980 


2.1195 


14.55 


3.90 


52.2 


1.0750 


2.0855 


12.30 


3.30 


80.3 


1.3420 


2.6365 


10.05 


2.69 


99.8 


1.2810 


2.5515 


7.80 


2.09 
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Table 2 continued 



r p (s) 


Xr (m) 




v\ (min 1 ) 


V2 (min 1 ) 


37.9 


0.6710 


2.0785 


14.55 


6.79 


46.9 


0.6640 


2.1185 


11.55 


5.39 


67.7 


0.7020 


2.1825 


8.55 


3.99 


108.4 


0.7360 


2.2665 


5.55 


2.59 


33.8 


0.3270 


2.1485 


14.54 


10.19 


45.3 


0.3360 


2.2295 


11.25 


7.88 


63.5 


0.3410 


2.2535 


8.10 


5.67 


103.3 


0.3340 


2.1745 


4.80 


3.36 



Table 3 

The results of the optimisation for fixed values of the parameter ji v from 0.05 to 
0.50s/m. The reduced \ 2 (x 2 /NDF) is also given (second column). The parameters 
A and B are correlated with (j, v ; of course, this is hardly surprising given the de- 
pendence of the tractive force on the final velocity of the vehicle (also remarked in 
Ref. [2]). A dependence of both A and B on the weight/load of/on the vehicle is 
expected (the results of the table correspond to a 16.6kg vehicle). The statistical 
uncertainties obtained from these fits are small when compared to the variation of 
the values with /jl v , at least in the interval of interest (0.10 < \x v < 0.20s/m). 



fi v (s/m) 


X 2 /NDF 


A (N-s/m) 


B(N) 


0.05 


1.414124 


1.987763 


4.836199 


0.10 


1.323325 


3.845995 


4.838278 


0.15 


1.302486 


5.706388 


4.840357 


0.20 


1.293971 


7.566957 


4.842438 


0.25 


1.289486 


9.427339 


4.844528 


0.30 


1.286761 


11.288190 


4.846608 


0.35 


1.284945 


13.148623 


4.848691 


0.40 


1.283655 


15.009084 


4.850781 


0.45 


1.282694 


16.868739 


4.852889 


0.50 


1.281952 


18.729476 


4.854964 
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Fig. 1. The mobile robot URANUS [5] utilising Mecanum wheels for omni-directional 
motion. 
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Fig. 2. Top view of a vehicle with four Mecanum wheels, along with its attached 
coordinate system (x,y). 
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Fig. 3. The driving (motor) force Fi acting on one of the wheels of the vehicle 
(chosen to be wheel 2), along with its decomposition into one force (Fi tP ) parallel 
to the rotational axis of the roller, which is in contact with the ground at that 
moment, and one in the transverse direction {Fi tj- The angle between the transverse 
component Fi t an d the rotational 'plane' of the wheel is denoted as a. The rollers 
shown are assumed to be those corresponding to the lower part of the wheel, part 
of which is in contact with the ground. 
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Fig. 4. Decomposition of the force of Fig. 3 into two forces F^ x and parallel 
to the axes of the attached coordinate system. The rollers shown are assumed to be 
those corresponding to the lower part of the wheel, part of which is in contact with 
the ground. 
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Fig. 5. Cancellation of the lateral forces in the special case uj\ = L02 = ^3 = W4. The 
rollers shown are assumed to be those corresponding to the lower part of the wheel, 
part of which is in contact with the ground. 
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Fig. 6. The attached coordinate system (x,y), along with the one having its origin 
at the centre of mass (CM) of the vehicle (x' ,y'). The rollers shown are assumed to 
be those corresponding to the lower part of the wheel, part of which is in contact 
with the ground. 
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Fig. 7. The relation between the laboratory coordinate system (X,Y) and the cen- 
tre-of-mass (CM) coordinate system (x',y') of the vehicle. 
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Position of the CM 
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Fig. 8. The position coordinates of the centre of mass (CM) of a Mecanum-wheeled 
vehicle in the laboratory coordinate system. The parameters generating this motion 
are given in Section 2.6. 
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Orientation angle vs time 




Fig. 9. The time dependence of the orientation angle (p. The parameters generating 
this motion are given in Section 2.6. 
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Position of the vehicle 




-60 -40 -20 20 40 60 80 100 

X (cm) 



Fig. 10. The position and orientation of a Mecanum-wheeled vehicle in the labora- 
tory coordinate system. The final position (i.e., at t = 4s) corresponds to the frame 
shown in black, the starting one in blue. The parameters generating this motion are 
given in Section 2.6. 
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Reproduction of the measured times 
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Fig. 11. Reproduction of the measured times needed for the accomplishment of 
motion for the 32 test cases of Table 2. The dotted line indicates equal measured 
(t ex P) and fitted (t th ) times. 
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